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1.  Introduction  and  Background  of  the  Project 

In  this  project  we  had  proposed  to  take  up  the  calculation  of  xenon  apparent  emission-excitation 
cross  sections  for  emission  lines  that  have  diagnostic  value  in  the  analysis  of  Xe-propelled 
electric  thruster  plasmas.  In  fact,  these  cross  sections  are  vital  for  models  that  are  necessary  for 
proper  integration  and  lifetime  predictions  of  Hall  effect  thrusters.  The  experimental 
measurement  of  cross  section  and  plasma  related  analysis  in  this  connection  have  been  in 
progress  and  are  being  carried  out  by  Dr.  R.  A.  Dressier  and  his  group  at  the  Air  Force  Research 
Laboratory,  Hanscom  AFB  [1],  To  supplement  this  work  we  had  planned  to  perform  at  I.I.T. 
Roorkee  in  India  reliable  calculations  for  electron  impact  excitation  cross  section  to  provide 
relevant  cross  section  data  needed  for  the  study  of  xenon  propelled  electric  thrusters.  Below  we 
briefly  introduce  the  xenon  propelled  electric  thrusters. 

An  ion  engine  relies  on  electrically  charged  atoms,  or  ions,  to  generate  thrust.  Xenon,  an 
inert,  noncombustible  gas,  is  electrically  charged  and  the  ions  are  accelerated  to  a  speed  of  about 
62,900  miles  per  hour  (30  kilometers  per  second).  The  ions  are  then  emitted  as  exhaust  from  the 
thruster  (see  figure  below  of  an  ion  engine  for  illustration),  creating  a  force,  which  propels  the 
spacecraft  in  the  opposite  direction.  The  primary  advantage  of  electric  propulsion  is  efficiency. 
An  ion  engine  is  10  times  more  efficient  than  its  alternative,  a  chemical  propulsion  system.  With 
xenon,  it  is  possible  to  reduce  propellant  mass  onboard  a  spacecraft  by  up  to  90  percent.  The 
advantages  of  having  less  onboard  propellant  include  a  lighter  spacecraft,  and,  since  launch  costs 
are  set  based  on  spacecraft  weight,  reduced  launch  cost. 


The  use  of  these  electrostatic  thrusters,  however,  poses  new  challenges  regarding  their 
proper  engineering  and  integration  on  satellites.  Precise  modeling  of  the  satellite  plasma 
environment  is  required.  A  suitable  approach  is  required  for  analyzing  xenon-propelled  Hall 
thruster  optical  radiation  based  on  apparent  electron  and  ion  impact  emission  cross  sections 
associated  with  lines  in  visible  and  near-infrared  region  of  the  spectrum. 

The  current  proposal  addresses  as  a  first  step  to  get  accurate  theoretical  electron  impact 
cross  sections  for  the  processes  where  initially  excited  metastable  states  viz.  I53  and  I55  of  xenon 
are  excited  by  electrons  to  the  2 pn  states  with  n  =1,10  in  Paschen’s  notation.  For  these  excitations 
no  other  theoretical  calculations  are  available  and  also  scarce  experimental  data  have  been 
reported  [2,3]. 

We  address  here  in  this  project  the  Hall  effect  thrusters  (HET),  which  are  a  high-specific 
impulse  alternative  to  chemical  propulsion  systems  of  spacecraft  [1,4].  In  HETs,  a  gas,  typically 
xenon  or  krypton,  is  efficiently  ionized  in  a  discharge,  and  positive  ions  are  electrostatically 
accelerated  to  generate  thrust.  The  performance,  plasma  characteristics,  and  durability  of  HETs 
are  extensively  evaluated  both  in  ground-based  test  facilities  and  in  space.  Basic  characteristics 
of  HET  plasmas,  such  as  charge  species  densities  and  temperatures  have  been  traditionally 
measured  by  various  plasma  probes.  However,  at  HET  conditions,  traditional  plasma-probe 
diagnostics  are  affected  by  problems  such  as  the  perturbation  of  the  local  environment  by  the 
probe,  the  complexity  associated  with  interpreting  the  probe  characteristics,  and  detrimental 
effects  in  regions  of  high  temperatures  (e.g.,  in  the  HET  discharge).  Consequently,  HET  plasma 
parameters  reported  from  probe  experiments  exhibit  significant  disparities.  Optical  plasma 
diagnostics  circumvent  these  problems  and  are  thus  an  attractive  alternative  to  probe 
measurements. 

A  key  component  of  collisional  radiative  model  (CRM)  is  a  set  of  emission  cross  sections 
associated  with  the  energy  transfer  processes  that  lead  to  population  of  states  that  radiatively 
relax  through  emission  in  the  spectral  region  of  interest  [1,4].  Proper  implementation  of  a  CRM 
has  been  primarily  hampered  by  the  lack  of  a  comprehensive  cross  section  set.  In  fact,  these 
cross  sections  are  vital  for  models  that  are  necessary  for  proper  integration  and  lifetime 
predictions  of  Hall  effect  thrusters  (HET).  If  we  focus  on  emission  excitation  cross  sections  for 
Xe  propelled  HETs,  the  optical  radiation  arises  from  electron  -  xenon  atom  collisions,  of  which 


the  most  important  contribution  can  be  assumed  to  be  coming  from  following  excitation 
processes: 

e_  +  Xe  — »•  Xe*  +  e_  (1) 

e"  +  Xem  -»•  Xe*  +  e"  (2) 

The  asterisks  in  processes  (1)  and  (2)  signify  excited  species  that  lead  to  optical 
emissions.  Superscript  m  indicates  metastable  states.  The  process  (1),  which  deals  with  the 
excitation  of  xenon  from  the  ground  state  to  upper  excited  states  have  been  studied  theoretically 
and  experimentally  [1-4].  The  process  (2)  which  deals  with  excitation  of  the  metastable  state  has 
not  been  explored. 

Metastable  states  of  the  noble  gases  are  important  constituents  of  plasmas.  In  addition  to 
their  long  lifetime,  these  states  have  large  inelastic  cross  sections  as  compared  to  the  excitation 
from  the  ground  state  [3].  Thus  it  is  important  to  include  these  processes  in  the  modeling  of 
plasma  which  is  prime  objective  of  the  present  project. 

2.  Outline  of  the  work  done 

•  Electron  impact  excitation  of  xenon  atoms  from  metastable  states  to  upper  2pn  (n  =  1-10) 
excited  states  have  been  considered. 

•  Excitation  cross  sections  are  calculated  using  the  Relativistic  distorted-wave  (RDW) 
theory. 

•  Fine-structure  atomic  states  are  represented  by  multi-  configuration  Dirac-Fock  wave 
functions. 

•  Continuum  projectile  electron  distorted  waves  are  obtained  through  Dirac  equations. 

•  Calculations  are  compared  with  experimental  results  [2]. 

•  The  cross  sections  are  incorporated  in  the  CRM  model  i.e.  in  the  analysis  of  Xe -propelled 
electric  thruster  plasmas,  being  carried  out  by  Dr.  R.  Dressier  at  AFB,  Hanscom  [1,4]. 


3.  Detailed  Description  of  the  Work  done 


3.1  Energy  levels  of  Xenon 

•  The  ground  state  of  xenon  has  5 p6  configuration. 

•  The  first  excited  state  has  a  5p56s  configuration  with  four  fine-structure  levels 
with  /=  0,  1,  1,  2. 

•  The  states  with  7  =  0  and  2  are  long-lived  metastable  states  i.e.  IS3  and  IS5  states. 

•  We  consider  the  electron  excitation  of  these  two  metastable  states  to  the  higher 
lying  5p56p  configuration  having  ten  fine- structure  levels  i.e.  2pn  (n  =  1-10) 
excited  states. 

The  figure  below  shows  the  energy  levels  of  the  xenon 
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Energy  levels  of  Xenon 

3.2  Description  of  the  states  in  relativistic  /-/  coupling  scheme 


The  ground  states  of  the  noble  gases  have  an  n p6  outer  shell  in  the  nonrelativistic  representation 
where  n  =  5  for  xenon. 

In  the  relativistic  j-j  coupling  scheme,  this  shell  is  broken  into  two  subshells  and  represented  as 
n  p  ~np  where  p  represents  a  p  electron  with  total  angular  momentum  (orbital  angular 
momentum  plus  spin)  j  =  1/2  while  p  has  j  =  3/2. 

Thus  the  metastable  states  have  configurations 

Paschen  notation 


5p25p36s  withJ=2 

5p5p46s  withJ=0 


ls5 

lS3 


Upper  states  have  configurations 
5p25p36p  with  J  =  0,  1,  2,  3 


5p25p36p  with/=  1,2 

5p5p%p  with  7=1,2 

5p5p46p  with  7  =  0,1 


These  upper  ten  states  are  represented  in  the  Paschen  notation  by  2pi  -  2p\o  in  order  of 
decreasing  energy  but  order  of  J  values  may  change  with  atom.  The  following  figure  shows 
diagrammatically  the  excitation  transitions  considered  in  the  present  work. 


Transitions  from  metastable  states  to  2pn  states 

of  Xe 


j 


4.  Calculation  of  Excitation  Cross  Sections  in  Xenon 

For  cross  section  calculation  we  use  relativistic  distorted  wave  (RDW)  method.  The  salient 
features  of  the  RDW  are  as  follows: 

4.1  Salient  features  of  RDW  method 

•  As  relativistic  effects  should  be  included  in  the  study  of  electron  impact  excitation  of 
heavy  atoms  where  spin-orbit  effects  are  important,  our  method  incorporates  full 
relativistic  effects. 

•  Our  method  considers  directly  the  fine- structure  target  states  without  any 
approximation. 


Our  method  includes  spin-polarized  electrons  in  a  natural  way. 


•  Our  Relativistic  Distorted  Wave  (RDW)  method  involves  solving  the  Dirac  equations 
to  describe  both  the  bound  and  continuum  electrons. 

•  The  relativistic  effects  are  thus  incorporated  to  all  orders  in  a  natural  way. 


4.2  Brief  Description  of  RDW  Theory 

The  distorted- wave  T-matrix  for  the  electron  impact  excitation  of  an  atom  having  N  electrons 
and  nuclear  charge  Z  from  an  initial  state  i  to  final  state  f  can  be  written  as  (atomic  units  are  used 
throughout) 

T™  =  {z~f  (1,  2, . ,N  +  l)\v~Uf(N  + 1)|  Aj,+  a  2, . .  N  + 1)) 

where  V  is  the  target-projectile  interaction  given  by 

Z  N  I 
rN+i  j=i  |rj  —  rN+1| 

here 

r j  (j  =  1,  . . .,  N)  position  co-ordinates  of  the  target  electrons 

Tn+i  position  co-ordinate  of  the  projectile  electron  with  respect  to  the  nucleus  of  the  atom 

Uf  distortion  potential 

A  antisymmetrization  operator 

The  wave  functions  X+1  ' ,  where  ‘ch’  refers  to  the  two  channels,  i.e.  initial  T  and  final  ‘f ,  are 
represented  as  a  product  of  the  N-electron  target  wave  functions  (jkh  and  a  projectile-electron 
distorted-wave  function  F°fy+(-)  ,i.e. 


XT  a,  2, . ,  N  +  l)  =  <(>ch(l,  2, . ,  N)  Fc°w+(_)  (kch ,  N  + 1). 


Here  ‘+’  refers  to  an  outgoing  wave  while  denotes  an  incoming  wave. 
The  projectile  electron  continuum  wavefunctions  are  given  by 


F±(k,r)= 


(2fl)3/2  Jfn 


±i'7, 


K  (k)— 
Km  y  >r 


'g^r)Z-an(f) 


E+c2 
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,1/2 
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1  //| 


j  m)Y*  ( k ) 
mi 


where 


(j.  magnetic  component  of  spin  of  projectile  electron 
kch  wave  vector  of  the  projectile  electron  in  the  initial  and  final  state 
f  larger  component  of  the  radial  part  of  the  wavefunction 
g  smaller  component  of  the  radial  part  of  the  wavefunction 
p  phase  shift 

E  relativistic  energy  of  the  incident  electron 


The  spin-orbit  functions  are  given  by 

W^>=£,  \vw 


and 

where 

J 

m 

V2)U<<T) 


X-Km<r,f)=  y  d  //'  V|  j  (r)i/i  UTI 


jLl  V 


£  ju 


total  angular  momentum  of  the  orbital 
the  z-component  of  j 

normalized  two-component  spin  wave  functions 


and  £=2  j-£ 

The  spin-angular  state  is  characterized  by  quantum  number  K,  which  is  defined  as 

if  j=t- 1/2  for  l  >0 

l-l  if  j=£+ 1/2 


K  - 


The  radial  parts  fK(r)  and  gK(r)  are  obtained  from  the  Dirac  equations 

i+fj  /d'-)T(c2-l,+v^r)-^V'i'’r)=0 

[ir-f}s^H^2-U+Ech>f^r)+±Wp(Kr)=0 

Here  U  is  distortion  potential  and  W  is  exchange  terms 
Dirac  equations  are  solved  using  the  boundary  conditions 

fK(r)- 
gK(r)- 


1  •  /,  I71  \ 
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Wave  function  On(h)  of  the  target  atom  is  solution  of  Dirac  equation  with  Hamiltonian 
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i= i 


-zcerV  +  /3c2  -- 
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N  i 

El2 


r  -r 
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The  single  electron  central  field  Dirac  orbitals  are  given  by 


O  (r  i 
r 


p„K(r) 

\  'QnJr)  X-k 


where 


a,  P  Dirac  matrices 

iVj  momentum  operator  of  the  jth  target  electron 

PnK  larger  component  of  the  radial  part  of  the  wavefunction 

QnK  smaller  component  of  the  radial  part  of  the  wavefunction 

The  bound  state  orbitals  satisfy  the  following  orthogonality  conditions, 

j  dr[PJ  r)PJ  r)  +  QJ  r)QJ  r)\  =  dn-n 
0 

(Xm  cr)  | &))  =  8k.k  Sm-m 


After  obtaining  the  wavefunctions  for  target  atom  bound  states  and  the  projectile  electron 
continuum  states.  We  evaluate  the  scattering  amplitude  for  the  excitation  of  a  noble  gas  atom  in  a 
metastable  state  with  total  angular  momentum  J,  and  magnetic  quantum  number  M\  to  a  final 
5p56p  state  with  angular  momentum  7f  ,  Mf  as 


f(Jf,Mf,VLf;Ji,Mi,VLi)  =  (2n)\U-T°wf(Jf,Mf,VLf;Ji,Mi,VLi) 


where  p,  and  pr  are  the  spin  projections  in  the  initial  and  final  channels.  Then  with  our 
normalization  the  differential  cross  section  (DCS)  is  given  by 


DCS=on,  ,  1 1  Zl f<JrMr\if;Jl,Ml,iLl  )\ 

£{  £J  i  -r  L  )  MiViMfPf 


and  the  ICS  is  obtained  by  integrating  the  DCS  over  all  scattering  angles. 


-! 


clQ. 


4.3  Calculation  of  the  Atomic  wave  functions 


The  metastable  states  are  single  configuration  wave  functions. 

<3>7=2  =  5p25p36s 

<£>y=0  =  5p5p46s 

The  excited  2 pn  states  are  linear  combinations  of  the  5p56p  configurations  with  the  same  J  value. 
=  cl5p25p36p  +  c25p25p36p+c35p5p46p +c45p5p46p 

The  mixing  coefficients  Cj’s  of  the  various  configurations  to  the  fine- structure  levels  of 
the  final  states  are  given  in  tables  1.  These  wave  functions  are  Dirac-Fock  wave  functions 
calculated  from  the  GRASP92.  To  optimize  the  wavefunctions  obtained,  we  have  compared  in 
table  2,  our  calculated  values  of  energy  differences  with  the  corresponding  experimental  values 
listed  in  NIST  Database  [5].  We  have  also  presented  in  table  3  the  comparison  of  our  calculated 
optical  oscillator  strengths  with  the  available  theoretical  [6]  and  experimental  [2]  results  for  the 
optically  allowed  transitions. 

Table  1.  Contributions  of  the  various  configurations  to  the  5p56p  states  of  Xe 


Configurations 

5  p  25p36  p 

5  p  25p36p 

5  p  5p46  p 

5  p  5p46p 

J  -0  levels 

2jPi 

0.0009 

0.9991 

2^5 

J  =  1  levels 

0.9991 

0.0009 

2p2 

0.0045 

0.0027 

0.0040 

0.9878 

2 Pa 

0.0001 

0.0019 

0.9937 

0.0043 

2 pi 

0.4306 

0.5679 

0.0015 

10‘7 

2pio 

J  =  2  levels 

0.5649 

0.4264 

0.0008 

0.0079 

2pi 

0.0006 

0.0004 

0.9990 

2 pe 

0.0527 

0.9467 

0.0006 

2p9 

0.9467 

0.0529 

0.0004 

Table  2.  Energy  differences  (eV)  for  the  transitions  considered:  NIST  -  experimental  values 
from  the  NIST  database  [5];  GRASP  -  calculated  values  using  the  GRASP92  program 


Transition  NIST  GRASP 


ls5-  2pl0 
ls5  -  2p9 
ls5  -  2pS 
lv5  -  2pl 
lv5  -  2p6 
ls5  -  2p5 
ls5  -  2p4 
ls5  -  2p3 
ls5  -  2p2 
ls5  -  2pl 

ls5-  ls3 


1.26483 

1.3703 

1.40542 

1.47398 

1.50577 

1.61816 

2.64229 

2.7394 

2.75383 

2.8259 

1.13187 


1.01583 

1.11376 

1.13561 

1.21215 

1.23393 

1.30729 

2.44226 

2.50259 

2.49999 

2.56711 

1.3239 


Table  3:  Dipole  oscillator  strengths  for  transitions  in  Xe:  GRASP92,  present  calculations; 
Theory,  theoretical  calculations  [6];  Expt.,  values  derived  from  high  energy  cross  section 
measurements  [2]. 


Initial  state 
Final  state 

GRASP 

ls5 

Theory  Expt 

1^3 

GRASP 

Theory 

2pio 

0.200 

0.268 

0.00021 

0.0249 

2p9 

0.112 

0.132 

2[h 

0.571 

0.641  0.59+0.19 

2 pi 

0.0227 

0.0117 

3. 09x1  O'5 

0.00496 

2p(, 

0.321 

0.253  0.30+0.10 

2p4 

6.09xl0'5 

7. 67x1  O'4 

0.482 

0.286 

2pi 

0.0009 

0.00266 

2pi 

0.004 

6.55xl0‘4 

0.760 

0.382 

5.  Xenon  Cross  Sections:  Results,  Discussion  and  Conclusion 


We  obsen’e  that  the  transitions  we  considered  are  of  the  following  type 


•  Transitions  from  the  \s$  J  =2  state  are  allowed  to  the  2 pn  states  with  J  =  1,  2,  3. 

•  Transitions  from  the  G3  /  =  0  state  are  allowed  only  to  the  2pn  states  with  7=1. 


•  All  other  transitions  are  forbidden  but  the  transition  from  I53  J  -  0  state  to  the  2 pn  states 
with  /  =  3  is  identically  zero  in  a  first-order  theory. 

The  figure  below  shows  diagrammatically  the  all  optically  allowed  transitions  from  metastable 
states  to  2pn  states  of  Xe 


Optically  allowed  transitions  from  metastable  states 

to  2pn  states  of  Xe 


J 


0 


y  2pr2p10 


2pu  y 


Fitting  of  the  calculated  allowed  transition  cross  sections'. 


For  allowed  transitions  we  find  the  cross  section  takes  on  the  Bethe-Bom  form  at  larger  energies 


G  -  \nal 


f 

EAE 


[In  (E)  +  b\ 


where 

E  energy  of  the  incident  electron  in  Rydbergs 
AE  excitation  energy  of  the  transition 
f  optical  oscillator  strength  of  the  transition 
b  fitting  constant 


We  have  fitted  our  excitation  cross  section  results  to  the  above  Bethe-Born  formula  and  the 
fitting  constant  ‘b’  for  the  allowed  transitions  is  given  in  table-4. 


Table  4:  Values  of  constant  b  for  allowed  transitions 


Atom 

Initial  State 
Final  State 

Xenon 

1^3 

1^5 

2p\o 

43.6192 

2.8511 

2p9 

- 

2.7067 

2 p& 

- 

2.6077 

2pi 

12.9306 

2.5839 

2p6 

- 

2.3348 

2 P5 

- 

- 

2 Pa 

2.4934 

1.3368 

2p3 

- 

1.2852 

2p2 

2.3252 

1.2952 

2/?i 

- 

- 

These  fitted  cross  sections,  thus,  allow  us  to  predict  the  cross  section  at  any  desired  energy  and 
can  be  easily  used  in  plasma  modeling  as  such. 

At  larger  energies  we  find  that  the  cross  sections  for  forbidden  transitions  behave  as  E'  where  E 
is  the  energy  of  the  incident  electron. 

Comparison  of  our  cross  sections  with  available  experimental  results 

Jung  et  al  [2]  have  recently  published  measurements  of  cross  sections  for  the  excitation  from  the 
metastable  I55  state  of  Xe  to  the  {2 p$  -  2pu})  manifold.  Only  two  of  these  transitions  ( 2p(,  and 
2 p$)  were  measured  at  higher  energies.  For  excitation  to  the  7=1  states,  which  are  allowed 
transitions  from  either  of  the  metastable  initial  states,  those  that  involve  a  change  of  core 
configurations  are  several  orders  of  magnitude  smaller  than  those  that  don’t.  We  note  that  the 
cross  section  for  the  IS5  -  2pj  transition  is  much  smaller  than  for  the  others  without  a  change  in 
core  which  is  in  conformity  with  the  magnitudes  of  the  oscillator  strengths. 

We  compare  our  results  with  the  measured  cross  sections  of  Jung  et  al  [2]  in 
Figures  1  and  2.  Our  cross  sections  agree  reasonably  well  with  the  measured  cross  sections  at 
larger  energies  for  the  excitation  of  the  2p^  and  2 p$  states  but  there  is  a  systematic  deviation  as 
the  incident  electron  energy  decreases  being  roughly  a  factor  of  two  larger  around  8  eV.  The 
agreement  with  those  transitions  which  were  not  measured  at  higher  energies  is  comparable  to 


these  two  cases.  Our  cross  sections  are  always  above  the  measured  values  at  8  eV  with  the 
exception  of  that  for  the  2 p5  state  which  may  imply  that  cascade  contributions  from  more  highly 
excited  states  are  larger  in  this  case  as  suggested  by  Jung  et  al  [2]. 


1  10  100 
Electron  energy  (eV) 


Figure  1:  Integrated  cross  sections  for  electron  excitation  of  xenon  from  the  I  .v5  state  to  the:  (a)  2 p6  state; 
(b)  2 p&  state.  Solid  curve,  presents  RDW  calculations;  experimental  points  with  error  bars,  Jung  et  al  [2], 


Figure  2:  Integrated  cross  sections  for  electron  excitation  of  xenon  from  the  I  .v5  state  to  the:  (a)  2 p5  state; 
(b)  2 p7  state;  (c)  2 pg  state;  (d)  2 pl0  state.  Solid  curve,  presents  RDW  calculations;  experimental  points 
with  error  bars,  Jung  et  al  [2]. 


The  following  conclusions  can  be  drawn  from  the  above  presentation  of  our  results 

•  The  RDW  method  has  been  shown  to  be  applicable  to  transitions  between  excited  states. 

•  Since  the  excitation  energy  of  these  transitions  is  relatively  small,  first-order  theories  are 
valid  at  lower  energies  than  for  excitations  from  the  ground  states. 

•  The  accuracy  of  the  cross  section  depends  crucially  on  the  accuracy  of  the  oscillator 
strengths  obtained  from  the  target  wave  functions. 

•  The  use  of  a  relativistic  formalism  (j-j  coupling)  clearly  explains  the  huge  variation  in  the 
magnitudes  of  these  cross  sections. 

•  We  have  presented  here  our  selected  results  and  these  are  now  published  [7]  and 
presented  at  various  conferences  (see  the  section  of  list  of  publications). 

•  Our  entire  results  published  and  unpublished  will  be  utilized  in  the  future  modeling 
calculations. 

•  We  have  achieved  our  target  and  obtained  the  aimed  metastable  excitation  cross  section 
results  of  xenon  required  for  the  CRM  model  of  Dressier  and  co  workers  [1,4]  and  these 
are  being  tested  and  new  plasma  modeling  results  are  under  calculation. 

6.  Application  of  Cross  sections  to  CRM  Model  of  Dressier  and  co  workers 

We  have  been  informed  that  our  calculations  of  electron-excitation  from  the  5p56s  7=2 
metastable  level  (lss  state  in  Paschen  notation)  to  the  lowest  six  5p56p  (2 p)  levels  have  allowed 
Dressier  and  co  workers  at  AFB  Hanscom  to  develop  a  collisional  radiative  model  (CRM)  for 
Xe  near-infrared  (NIR)  emissions  based  on  population  and  depopulation  of  the  metastable  level 
through  I55  -  2 p  transitions.  Application  of  the  model  to  spectral  intensities  observed  in  the 
plume  of  a  Hall  thruster,  however,  demonstrates  that  the  model  overestimates  the  population  of 
the  1  s 5  state.  The  observed  spectrum  can  be  reproduced  with  high  fidelity  if  additional 
metastable  de-excitation  mechanisms  are  effective. 

On  the  demand  of  the  Hanscom  Group,  we  have  further  carried  out  calculations  for  the  excitation 
of  xenon  from  its  I55  level  to  1  s%,  l.s’3  and  l.v4  levels.  Thereafter,  the  group  has  demonstrated  that 
the  electron-induced  depopulation  through  the  excitation  to  5p56s  7=1  ( 1  ,v4)  and  other  5p56s 
levels,  for  which  newly  calculated  cross  sections  are  presented,  can  account  for  the  additional  de- 


excitation  mechanisms  and  give  reliable  results.  This  work  we  are  planning  to  present  in  the 
GEC-2007  meeting  to  be  held  at  Arlington  VA,  during  October  2  -5,  2007  (see  list  of 
publications). 


7.  Future  Work 

In  the  light  of  recent  experimental  results  of  Jung  et  al  [8]  we  have  carried  out  relativistic 
distorted  wave  calculations  for  electron  impact  excitation  of  the  ten  higher-lying  fine- structure 
levels  of  the  3p55p  configuration  of  argon  from  the  lowest  metastable  states  (the  J  =  0,2  levels  of 
the  3p5  4s  configuration).  We  compare  our  theoretical  results  with  their  experimental  results  and 
discuss  the  differences  from  the  similar  excitation  to  the  3p54p  levels  from  the  same  metastable 
states  [9].  In  the  light  of  this  work  on  argon,  we  feel  that  in  xenon  also  the  electron  impact 
excitation  of  the  ten  higher-lying  fine-structure  levels  of  the  5ps7p  configuration  from  the  lowest 
metastable  states  (the  J  =  0,2  levels  of  the  5p56s  configuration)  possibly  can  contribute.  This  has 
to  be  tested  and  we  are  currently  working  on  this  aspect  as  well.  In  addition,  as  pointed  out  by 
the  Hanscom  group  that  the  electron  excitation  of  Xe+  ion  can  influence  the  results  of  CRM 
model,  we  would  therefore,  like  to  test  this  as  well.  Currently,  we  are  also  planning  to  extend  our 
RDW  calculations  to  the  excitation  of  Xe+  ions  which  will  be  our  future  work. 
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